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I examine the role of final state interactions in cold nuclear matter in modifying hadron production 
on nuclear targets with leptonic or hadronic beams. I demonstrate the extent to which available 
experimental data in electron-nucleus collisions can give direct information on final state effects in 
hadron-nucleus and nucleus-nucleus collisions. For hadron-nucleus collisions, a theoretical estimate 
based on a parton energy loss model tested in lepton-nucleus collisions shows a large effect on 
mid-rapidity hadrons at fixed target experiments. At RfflC energy, the effect is large for negative 
rapidity hadrons, but mild at midrapidity. This final state cold hadron quenching needs to be taken 
into account in jet tomographic analysis of the medium created in nucleus- nucleus collisions. 



I. INTRODUCTION 

Hadron production on nuclear targets is strongly in- 
fluenced by the presence of cold and hot nuclear matter. 
The most spectacular effects is jet quenching in nucleus- 
nucleus [A + A) collisions [1-4], namely, the suppres- 
sion of hadron production at large transverse momentum 
compared to a suitably scaled cross section in proton- 
proton collisions. This phenomenon is widely used as a 
tool to explore the properties of the hot QCD medium 
produced in heavy-ion collisions [5-7], and as evidence 
for the creation of a novel state of matter at the Rel- 
ativistic Heavy Ion Collider (RHIC) [1-4], possibly the 
Quark-Gluon Plasma (QGP) [8]. Hadron suppression has 
also been observed in lepton-nucleus {^ + A) [9-14] and 
hadron-nucleus {h + A) collisions [15-17], where one does 
not expect the formation of an extended hot medium. In 
this case, the target nucleus itself (cold nuclear matter) 
induces the observed suppression of hadron production, 
which I will refer to as "cold hadron quenching" , or cold 
quenching in short. 

Nuclear effects in cold nuclear matter can be classi- 
fied as initial-state or final-state depending on whether 
they happen before or after the hard collision which gen- 
erates the hard probe. Final state (FS) effects can be 
isolated in semi- inclusive hadron production in i + A col- 
lisions. Hadron suppression in these collisions [9-13] is 
typically attributed to radiative energy loss of the struck 
quark or to nuclear absorption of a colorless prehadron, 
see [18] and references therein. Initial state (IS) effects 
can be experimentally isolated in Drell-Yan processes in 
h+A collisions [19-21], and they are attributed to nuclear 
shadowing or radiative energy loss of the incoming parton 
[22-24]. In large-pT hadron production vah + A collisions 
both IS and FS effects are present, and they cannot be 
easily disentangled. They give rise to a host of interest- 
ing effects. They modify the shape of midrapidity hadron 
Pt spectra, slightly suppressing it at small pT I — 2 
GeV and enhancing it at intermediate 2 GeV < Pt ^ 6 
GeV; this is known as Cronin effect [25-27]. An exten- 
sive study of the rapidity dependence of the Cronin effect 
in deuteron-gold {d -t- Au) collisions at RHIC has shown 
an interesting evolution of the Cronin effect: the Cronin 
enhancement grows at backward rapidity y — ycm < 



[28-31]; however, at forward rapidity y — ycm > 0, the 
spectrum is suppressed up to large px [15-17], a trend 
confirmed also at lower energy collisions [32, 33]. The 
forward rapidity suppression has been widely interpreted 
as evidence for the onset of the Color Glass Condensate, 
a universal state of highly saturated quark and gluons in 
the nuclear wave function [34, 35]. However, explanations 
in terms of IS energy loss and higher-twist shadowing 
[33], leading- twist shadowing [36], Sudakov suppression 
[37] or FS parton recombination [38] have been equally 
successful in describing the data. The rise of the Cronin 
effect at backward rapidity is more difficult to under- 
stand [39]. Explanations in terms of IS anti-shadowing 
[40] or saturation [28] have been proposed. Finally, re- 
cent PHENIX data on neutral pion production in c? -I- Au 
collisions at midrapidity suggest a small suppression of 
midrapidty tt" at px > 10 GeV [41]. It cannot be ex- 
plained by the EMC effect, which is effective at ^ 15 
GeV [42, 43], but may accommodate a small final state 
energy loss of order 10% [43]. 

A consistent interpretation of this wealth of experi- 
mental data requires a deep understanding of IS and FS 
interactions at the parton and hadron level, and the de- 
velopment of a unified computational framework [44] . As 
a contribution to this program, in this paper I will an- 
alyze phenomenologically the contribution of final state 
interactions to hadron production in h + A and A + A 
collisions, and I will show that it is indeed non negligible 
in the whole backward rapidity hemisphere up to RHIC 
energy. At the Large Hadron Collider (LHC) it will be 
important only at very backward rapidity y^ycm ^ 3. In 
Section II, I will review the kinematics of hadron produc- 
tion at leading order (LO) in perturbative QCD for Deep 
Inelastic Scattering (DIS) and nucleon-nucleon collisions 
(NN), and build a dictionary that relates the kinematic 
variables used in the 2 cases. In Section III, I will com- 
pare the NN and DIS phase spaces at present and future 
experiments in terms of either set of variables, to under- 
stand in detail the relevance of NN to DIS and vicev- 
ersa. In Section IV, I will use the developed kinematic 
dictionary to show the extent to which present i -\- A ex- 
perimental data can give information on final state cold 
nuclear matter effects in h + A and A + A collisions (col- 
lectively, A + B collisions) . Then, I will use an energy 
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FIG. 1: Left: LO kinematics for parton production in NN 
collisions - double lines indicate hadrons or nuclei, single lines 
are partons. Right: LO kinematics for parton production in 
DIS collisions - double lines indicate hadrons or nuclei, thin 
single lines are partons (bottom) or leptons (top). The labels 
define the particles 4-momenta. 



loss model tuned to £+A data to estimate cold quenching 
in h + A collisions for midrapidity hadrons at the Super 
Proton Synchrotron (SPS) and at Fermilab (FNAL), and 
at various rapidities at RHIC. A comparison of PS and IS 
effects will be discussed in Section V, and my conclusions 
reported in Section VI. 



II. PARTON PRODUCTION IN DIS AND NN 
COLLISIONS 

Considering parton and hadron production at LO in 
NN collisions and DIS collisions, it is easy to provide an 
explicit dictionary translating between the variables tra- 
ditionally used in the analysis of the 2 processes. I will 
start by discussing the kinematics of parton and hadron 
production in NN collisions in the center-of-mass frame 
(c.m.f.). I will then suitably identify the DIS kinematics 
in terms of NN variables, and derive the dictionary. I will 
work in the framework of collinear factorization in pQCD, 
and use light-cone coordinates throughout this discus- 
sion: for any 4- vector I write — (a"*", a~, ot), where 
{a^ ±a^)/y/2 are the plus- and minus-momenta and 



flT = the transverse momentum. 



NN collisions 



In pQCD at leading order in the coupling constant 
as , parton production in NN collisions proceeds through 
2^2 partonic collisions (see Fig. 1 left and Table II for 



the definition of kinematic variables.) Several LO pro- 
cesses can contribute to a given ij—^fif2 collisions, rep- 
resented by a black disk in the cartoon, see Ref [45] for 
details. The momenta of the 2 nucleons colliding in the 
c.m.f. with energy y/s/2 each are 



/ = 



J = 



2 

V2S' 



/2S 



,0t 



(1) 



where M is the nucleon mass and 



1 + v/l + MVs^ 



(2) 



I will neglect terms of order 0{M'^ /s) compared to terms 
of 0(1), and will use s « s. Note also that in the defini- 
tion of the nucleon momenta, I explicitly retain the nu- 
cleon mass in Eq. (1) to be able to perform boosts to the 
rest frame of either nucleon. If we assume the partons to 
be massless and collinear to their parent nucleons, their 
4-momenta in terms of the parton fractional momenta Xi 
read 



I = (a;i — ,0,0t 



J 



(0,X2^,0t 



(3) 



In terms of rapidities and transverse momentum pT, the 
parton 4-momenta read 



(4) 
(5) 



We can express the parton fractional momenta in terms 
of PT, Vi as 



X2 = ^(e-^i +e-y^) 



(6) 



Finally, the Mandelstam invariants are defined as follows, 

s = (z + if 



i={^-h? = {h-3f 
n={i-h? = {h-3f 



(7) 



and 4-momentum conservation is expressed as J + i + m = 
0. In terms of rapidities and transverse momentum, the 
Mandelstam invariants read 



S — X\X2S 
„2 



(8) 
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Variable Defimtion 
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Xl 


= i+/7+ 


Initial state projectile parton fractional momentum. 


X2 


— 01'' 


lllililcLi otcLtc tcLlgcb pdl toll il dC tlOlldil lllOlllclltLllll. 


PiT 


= \fiT\ 


Final state partons transverse momentum (relative to beam). 


Vi 


= 0.51og(/+//-) 
= 0.5 log (1^) 


Final state partons rapidity. 


ycm 


Rapidity of the center of mass. 


z 


= Pt/f^ 


Hadron fractional momentum relative to parent parton /i. 


PhT 


= \PhT\ 


Hadron transverse momentum (relative to beam). 


Vh 


= 0.51og(p+/p-) 


Hadron rapidity. 


n 


= -logtan(r/2) 


Hadron pseudorapidity {6* is the angle between the parton 
and the beam in the center of mass reference frame) . 



TABLE I: Definitions of the kinematic variables for semi-inclusive parton and hadron production in pQCD (top and bottom 
part of the table, respectively). Particle 4-momenta are defined in Fig. 1. Boldface variables are experimentally measurable. 
The remaining variables are theoretically defined in the QCD parton model in coUinear factorization. Note that at LO, with 2 
final state partons, pit = —P2T = Pt- 



Variable Definition Target rest frame 





= P2 




Target mass. 


Xb 


_ 

2P-q 

= -q^ 




Bjorken scaling variable. 






Negative four-momentum squared of the virtual photon. 


V 


— qP 




Energy of the virtual photon in the target rest frame. 


V 


_ qP 

k-P 


1/ 


Fractional energy loss of the incident lepton. 




= (P + lf 




Invariant mass squared of the hadronic final state. 




_ Ph-P 

qP 


= Sh 


Fraction of the virtual photon energy carried by the hadron. 


Pt 


= \Pt\ 




Hadron transverse momentum (relative to the virtual photon momentum) . 



TABLE II: Definitions of the kinematic variables for semi-inclusive DIS. The Lorentz invariant definition and its form in the 
target rest frame are provided. Particle 4-momenta axe defined in Fig. 1. All variables are experimentally measurable, hence 
typeset in boldface. Note that xb = /{2Mv) independently of the chosen reference frame. 



In order to compare collider and fixed target experiments, 
and different beam energies, it is useful to consider the 
rapidity in the c.m.f.: 

yc.m.f. y Vcm ■ (9) 

The backward rapidity region (target hemisphere) corre- 
sponds to y — ycm < 0, and the forward rapidity region 
(projectile hemisphere) to y — ycm > 0. 

Hadronization in the collinear factorization framework 
proceeds through independent parton fragmentation into 
a hadron. It is universal, i.e., independent of the process 
which produced the fragmenting hadron, e.g., NN or DIS 
collisions [46]. The hadron fractional momentum z is 
defined by 

Ph = zfi 

^ i (10) 

PhT = zfiT . 

Therefore the on-shell hadron momentum ph reads 

Ph^{zf^,'^^^^^^,zhT). (11) 

The parton and hadron rapidities are related by yi = 
yh + iog{mhT / Phx) ■ The non perturbative dynamics of 
the fragmentation process is encoded in universal frag- 
mentation functions, which can be obtained in global fits 
of experimental data [47, 48] . 



B. DIS collisions 

At LO in pQCD, deeply inelastic scattering proceeds 
by exchange of a virtual photon in the f-channel, explic- 
itly shown in Fig. 1 right. The DIS Lorentz invariants are 
defined in Table II. Semi-inclusive nDIS is best discussed 
in terms of u and Q^, which are the most relevant vari- 
ables to hadron quenching processes in nuclear targets. 
Analysis of inclusive DIS is usually carried out using x b 
and Q'^. 

DIS experiments can be performed with a fixed tar- 
get (ft) or in collider mode (cl). Examples are the EMC, 
HERMES, JLAB experiments, and the Electron-Ion Col- 
lider (EIC), respectively. The colliding nucleon and lep- 
ton momenta are 

Pft = ^:0t) , fc/t = (x/2£;e,0,0T) 

Pel = , V2En , Ot) , kci = (y2Ec , 0, Ot) 

(12) 

where and Ej^ are the electron and nucleon energies 
measured in the laboratory frame. To discuss both modes 
at the same time, it is convenient to introduce the target 
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rest frame energy of the electron, El^^: 



E, 



Ef=[2LE. (13) 

Then the invariant y for both modes becomes y = v/E^L 



fixed target 
coUider mode 



C. The dictionary 

To connect the DIS and NN kinematics, wc can boost 
the DIS collision to a frame in which the target has en- 
ergy \fs/2 per nucleon. Then, we can imagine the lepton 
to be a parton of a phantom miclcon of energy \fs/2 
collinear with the lepton, with 4- momentum P'^ = P^. 
Comparing the top and bottom of Fig.l we can identify 



P = J, P' = 7, k. 



k' = f2 



(14) 



The virtual photon momentum q, the fractional momen- 
tum Xe of the initial state lepton and the rapidity ye of 
the final state lepton arc identified as follows 



k-k' 



/2 



= fc+/P' 



Xl, ye 



2/2 



In this way, we can relate the DIS kinematics to the NN 
kinematics discussed in Sect. II A. As an example, it is 
immediate to sec that, in terms of NN variables, = —t. 
The full translation dictionary from DIS to NN variables 
can be obtained in a straightforward way by combining 
the results of Sects. IIA-IIB and the definitions of Ta- 
bles I-II. 

First, we can express the DIS invariants in terms of 
parton rapidities and transverse momenta. Neglecting 
target mass corrections, i.e., up to terms of 0{M^/s), we 
obtain 

XB = ^{e-y^ + e-y^) 



2M 



(15) 



1 



y = 



1 -I- e?/2-yi 
Zh = z . 

Note that the first 3 variables are not independent be- 
cause = 2Mxbi^, and that xb = X2 is interpreted as 
the struck parton fractional momentum, as expected in 
DIS at LO. Note also that i/ increases with increasing px 
and increasing yi. In other words, a parton of positive 
and large yi travels in the opposite direction as its par- 
ent nucleon, hence in the target rest frame it is very fast. 
Conversely, a parton of negative and large yi travels in 
the same direction as its parent nucleon, which means 
quite slow in the target rest frame. It is also interesting 
to note that up to terms of order 0(M^/s), the parton 
and hadron energy in the target rest frame are 



£ = 1/ Eh = ZhV . 



(16) 



Finally, we can invert Eq. (15) to obtain the NN variables 
in terms of DIS invariants: 

= (1 - y)Q^ 



yi 



= - log ( 



2MEf^ y 
l-y^ 



2/2 = 2/1+ log (~^) 

Z = Zh 



(17) 



with y = v/Ep. 

Note that in DIS, the electron energy E^^^, hence the 
electron Xe, is fixed by the experimental conditions: this 
is different from NN collisions where the parton j has 
an unconstrained fractional momentum. Changing the 
c.m.f. energy to \/s' simply results in a shift of the parton 
rapidity, 



2/1 



2/1 + Ayi 



(18) 



where Ayi = \og{y/s/Vs'). The value of Ayi compared 
to RHIC top energy ^/s = 200 GeV is listed in Table III 

for the experiments of interest in this paper. Another 
difference between DIS and NN collisions is the rapid- 
ity difference Ay between the outgoing "partons" . In 
DIS, the electron fractional momentum is fixed, so that 
At/|£)75 = yi^Ve = log {y/{l~y)) is determined for each 
Pt and yi by the corresponding value oi y = v/E^J'^, and 
can span only a limited range: 



1 / ymin \ ^ A ^1 f Vmax 

log -; < Ay\Dis < log 



(19) 



For example, at HERMES the experimental acceptance 
0.07 < y < 0.85 translates into -1.1 < Ay^Dis < 0.75. 
In NN collisions, neither parton fractional momentum 
is fixed by the experimental conditions, hence Aj/jviv = 

2/1 — 2/2 can span 



■ log (- 



Pt 



l) < Ay\NN < log ( 



y/se 



Pt 



(20) 



For example, for an observed parton with yi = —2 and 
Pt = 2 GeV, corresponding to the middle of the HER- 
MES DIS phase space, we obtain —2.5 < Ay|jvjv < 6.6 , 
even though the average (Ay|jvjv) will lay in the middle 
of this interval. 



III. COMPARING THE PHASE SPACES 

We can now compare in detail the phase spaces for 
parton production in NN and DIS collisions. For this 
purpose, I will define a NN-equivalent DIS phase space 
and a DIS-equivalent NN phase space. 
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FIG. 2: Left: the RHIC-equivalent phase space of DIS experiments at = 27.6 GeV (HERMES, sohd hne), at = 12 GeV 
(HERMES and JLAB, dashed hne), and at Ee = 280 GeV (EMC, dot-dashed line). The dotted line shows the borders of 
the LO pQCD phase space at top RHIC energy, -^s = 200 GeV. The 2 arrows show the location of the midrapidity region 
at SPS and FNAL fixed target experiments. The open and filled symbols show the position of a representative selection of 
data on hadron suppression collected at the HERMES [13, 49] and EMC experiments [10], respectively. Red squares are for 
Zh-distributions, green triangles for (/-distributions and blue circles for Q'^-distributions. The systematic uncertainties on EMC 
data, shown as error bars, are discussed in Section IV A. Right: NN-equivalent EMC and COMPASS phase space at -^s — 27.4 
GeV, compared to the SPS and FNAL phase spaces. 





SPS FNAL RHIC RHIC LHC 


Vs [GeV] 
Ayi 


17.5 27.4 63 200 5500 
2.4 2.0 1.2 -3.3 



TABLE HI: Rapidity shifts Ayi of the RHIC-equivalent DIS 
phase space, tabulated for some energies of interest. 



and y„ 



■' min ' min 



Q 



2 

min 



< 



2M 



2Mv - Wi 



(21) 



A. NN-equivalent DIS phase space 



Given a DIS phase space, i.e., a given experiment ac- 
ceptance region in the {v, Q^) plane, I define its NN- 
equivalent phase space as its image in the {pT,yi) under 
Eqs. (17). (I do not consider the transformation of the 
fragmentation variable Zh into z because of its triviality.) 
The reason for this definition is that for both NN and 
DIS collisions we can identify the parton /i of Fig. 1 with 
the "observed" parton in NN and DIS collisions, i.e., the 
parton which fragments into the observed hadron. Then 
the variables pr and yi fully characterize the observed 
parton. An analogous definition holds when using xb 
instead of u as independent variable. 

As an example, the HERMES DIS phase space in the 
(i/, Q^) plane is determined by the values of W^j„, (3^i„ 



Additionally, one may impose stronger cuts on i>, e.g., 
> i^rnin, as at the EMC experiment, and in some HER- 
MES analysis. 

With Eqs. (17) it is easy to plot the NN-equivalent 
DIS phase space in the (yi,PT) plane. As an example, 
we can consider the RHIC-equivalent phase space of the 
HERMES and EMC experiments, using = 200 GeV, 
shown in Fig. 2 left. Note that according to Eq. (18), 
the NN-equivalent phase space at other center of mass 
energies can be obtained by a shift yi^yi + Aj/i, see 
Table III. I assume the pQCD formulae used to define 
the NN-equivalent phase space to be valid at RHIC top 
energy for p > po = 1 GeV: the corresponding pQCD 
confidence region is plotted as a dotted line, see Eq. (23) 
for details. 

We can see that the HERMES experiment, with^E*'"-'^ = 
12 and 27.6 GeV, covers less than one third of the avail- 
able RHIC pt range at yi ~ —3, with shrinking px cov- 
erage at larger rapidity. In the SPS/FNAL midrapidity 
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region it reaches pt = 2.5 GeV at most. Since 



yi < log ( 



2MET^ymax' 



(22) 



and Umax cannot be increased above 1, the only way to 
effectively reach larger values of yi is to increase the elec- 
tron beam energy E^^^. Indeed, the EMC experiment, 
with E^^f = 100 — 280 GeV, covers a larger span in ra- 
pidity and extends to yi > 0. Moreover, the increased 
energy allows in principle to reach much higher px than 
at HERMES. However, only the < 3 GeV region has 
been explored. As also shown in Fig. 2 left, the proposed 
Electron-Ion Collider (EIC) [50, 51] will be able to ef- 
fectively study the yi > region, and cover most of the 
RHIC phase space. Likewise, it will cover only the yi < 
part of the LHC phase space. 

The reason why present experimental data ini+A col- 
lisions reach only px ^ 2 GeV, is that conventional DIS 
variables Zh, v or explore the available NN-equivalent 
phase space in an uneven way. Moreover, in single differ- 
ential distributions like dN^/dz, the integration over the 
remaining variables favors low values of Q^, hence low- 
Pt values. While HERMES is inherently limited in its 
Pt reach by the low electron beam energy, the EMC ex- 
periment covers, in principle, most of the SPS and FNAL 
phase space, see Fig. 2 right. Therefore, a rebinning of 
the EMC experimental data in terms of NN variables 
would result in an experimental measurements of final 
state nuclear effects, much needed for correctly interpret- 
ing large-pT hadron spectra in h + A and A -I- yl collisions 
at SPS and FNAL. Another possibility would be to study 
fi + A collisions at the COMPASS experiment [52], which 
has a muon beam energy of E^ = 160 GeV comparable 
to EMC, and whose phase space is also shown in the plot. 



B. DIS-equivalent NN phase space 

When discussing NN collisions in the framework of 

coUinear factorization in pQCD, we should first define the 
region of validity of perturbative computations: pt > Po- 
Typically one needs ^ 1 GeV, which agrees with the 
phenomenological analysis of Refs. [42, 53]. Then, the 
NN phase space at a given yi is defined by the kinematic 
bounds on 2— >2 parton scatterings [42]: 



|yi| < cosh 
Po<Pt< 



y^po) 



2cosh(2/i) 

- log - e-^^) < t/2 < log - e^^) 

\PT ' ^PT ' 



(23) 



< z<l 



Introduction of intrinsic parton transverse momentum in 
the formalism, or use of next-to-leading order kinematics 



[54] , would relax somewhat these bounds. We should also 
keep in mind that at large rapidity, where the 2— >2 phase 
space is becoming more and more restricted, 2— >1 parton 
fusion processes may become the dominant mechanism, 
because they are sensitive to much lower fractional mo- 
menta Xi [39]. Hence, at the boundary of the NN phase 
space, the presented analysis becomes unreliable. 

The DIS-equivalent NN phase space is defined as the 
image of Eqs. (23) in the (y,Q'^ ,y, Zh) space under 
Eqs. (15). It is 4-dimensional and difficult to directly vi- 
sualize. A way around this problem is to define suitable 
trajectories in NN phase space averaged over 2/2, and to 
project them into the DIS-equivalent {i',Q^) and {i^,Zh) 
phase spaces. We can define a PhT- and y/i-dependent 
average observable as follows 



(O) 



J dz dyi dy2 0{pT,yuy2, z)-^ 



dp^dy\dy^dz 



PhT,Vh 



Jdzdyidy2-S 



a^AB^hX 



where 



da 



AB^hX 



dp'^dyidy2dz 



E 

/i 



1 . do-^^^/i^ 

-2D)^{^) \, . 
z'^ " dp^dyidy2 



(24) 



(25) 



^^AB^fiX ^Yic LO pQCD differential cross-section for 
production of a /i parton pair in a collision of hadrons A 
and B (nucleons or nuclei), and is its fragmentation 
function into the observed hadron, see Ref [42] for details 
[80]. Then, we can use Eqs. (15) to compute {v)pf^T,yh^ 
{Q'^)phT,yhJ and {zh)phT,Vh- These values parametrize the 
fixed-y/i trajectories {{{v)pT,y, {Q'^)pT,y)\PT > Po] and 
{{{v)pT,yi {zh)pT.y):PT > Po} in the DIS-equivalent phase 
space. 

As an example, in Fig. 3 I considered NN collisions 
at RHIC top energy ^/s = 200 GeV and at fixed tar- 
get energies y/s = 17 — 27 GeV, and plotted the fixed- 
yh trajectories in the DIS-equivalent phase space. The 
range of pt spanned along each trajectory is tabulated 
in Table IV. The spanned range in is limited by the 
maximum pT at each rapidity, according to Eq. (23). As 
expected, the larger the rapidity yh ~ yi the smaller 
the spanned u. RHIC trajectories with yh ^ —'^ span 
pretty low values of 1^ < 60 and large values of Zh > 0.5, 
where the EMC and HERMES experiments have shown 
non negligible cold nuclear matter suppression of hadron 
production. At higher rapidity, the larger spanned values 
of v will make cold nuclear matter effects less prominent. 
The consequences of these remarks for the interpretation 
of hadron production in h + A and A + A collisions will 
be further discussed in Section IV and V. 



IV. FINAL STATE COLD NUCLEAR 
QUENCHING IN h -|- A COLLISIONS. 

As we have seen, a parton produced at negative ra- 
pidity, y — ycm < 0, in a, h + A collision travels in the 
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FIG. 3: Right: Fixed-yi NN trajectories plotted in the DIS-equivalent {v,Q^) phase space for RHIC at = 200 GeV and 
various rapidities, for FNAL and SPS at midrapidity. The dot-dashed line encloses the HERMES phase space; the dashed line 
encloses the EMC phase space. The arrow indicates the direction of increasing (pr) and (zh). Left: Trajectories in the (v,Zh) 
plane. The arrows indicate increasing pr and Q2. 





SPS 


FNAL 


RHIC 




= 17.5 GcV 


= 27.4 GeV 


^ = 200 GeV 


Vh 








-1 -2 -3 


PhT 


1-8 


1-12 


1-90 1-60 1-25 1-9 



TABLE IV: Range of average {pT) spanned along RHIC tra- 
jectories at fixed rapidity y\ and ^/a — 200 GeV. phT is quoted 
in GeV. 



hard 

IS parton int. FS partons FS prehaclron hadrc 



3, .9 




same direction as the target nucleus: seen in the nucleus 
rest frame, it appears to move slowly and corresponds 
to a low value of ly in the language of ^ -I- ^ collisions. 
Therefore, based on the observed suppression of hadron 
production in lepton-nucleus DIS [9-12, 12] at low i^, and 
on the kinematic analogy between DIS and NN collisions 
discussed in the previous sections, we can expect non- 
negligible hadron suppression due to FS interactions in 
cold nuclear matter also in h + A and A + A collisions. 

Discussion of medium effects is best carried out in the 
medium rest frame: in the case of cold nuclear matter in 
£ + A and h + A collisions it is the nucleus rest frame. 
I am interested here in processes characterized by large 
values of xb = 2^2 ^ 0.1, typical of semi-inclusive nDIS 
measurements at HERMES and large px hadron produc- 
tion at not too forward rapidity in h + A collisions. In 
this regime, the hard interaction is well localized inside 
the nucleus, and the nucleons act incoherently as targets 
[55]. The process evolves in time as follows, see Fig. 4. 
First the electron (or a parton belonging to the proton) 
penetrates the nucleus, and undergoes a localized hard 
collisions. Then, a "final-state" system of 1 electron and 
1 parton (or 2 partons) is produced, with both particles 
essentially traveling along the beam direction, even for 
rapidity values far from the center of mass rapidity in the 
target hemisphere {y— ycm < 0). Later on the final state 



hard 

IS lepton int. FS jjaitoirs FS preliadron hadroir 




FIG. 4: Top: Initial and final state interactions in h + A 
collisions in the nucleus rest frame. Bottom: absence of initial 
state interactions in £ + A collisions. The nucleus is drawn as 
an oblong oval for convenience only. 



partons hadronize and one of the produced hadrons is 
detected. The time scale on which hadronization starts 
after the hard interaction is not at present well known 
[18]; it may be as small as the nuclear radius [13, 56, 57], 
in which case the hadronization process would start in 
the medium. Nuclear medium effects may be classified 
as initial state (IS) effects on particles before the hard in- 
teractions, or final state (FS) effects on particles created 
after the hard interaction. In the case of £ + A colli- 
sions, electromagnetic reinteractions of the incoming or 
outgoing lepton are suppressed compared to the strong 
FS reinteraction of the parton and hadronizing system. 
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For h + A collisions one needs in principle to account for 
both IS and FS interactions [44]. 

Except at very forward rapidity, I will assume FS and 
IS effects in h + A collisions to be factorizable because 
of the large rapidity difference between the IS and FS 
partons induced by the hard scattering. Then, I will ex- 
plore the possible size of FS effects on single inclusive 
hadron production. Differently from £ + A collisions, the 
proton projectile interacts with the nucleons along its 
trajectory. The hard parton produced in the hard col- 
lision starts propagating at nearly the speed of light in 
the same direction but slightly behind the projectile pro- 
ton. The time scale for parton production, thard oc 1/Q 
is miich smaller than the time scale for soft particle pro- 
duction in proton-nucleon collisions, tsoft oc 1/Aqcd- 
Hence, we may assume the nuclear medium traversed by 
the produced parton in h + A collisions to have approx- 
imately the same properties as the cold nuclear matter 
which would be traversed in i + A collisions, i.e., the tar- 
get nucleus itself. With this in mind, we may assume 
final state hadron quenching effects to be comparable in 
the 2 cases. 



A. Lessons from £-\- A data 

Nuclear modifications of hadron production in £ + A 
collisions arc typically studied in terms of the hadron 
multiplicity ratio 



RM{^h,nu,Q'^) = 



1 



1 



dzhdvdQ"^ 



iV^^^ dzhdvdQ"^ ' 
(26) 



i.e., the single hadron multiplicity on a target of mass 
number A normalized to the multiplicity on a deuteron 

target. Then, we can use the dictionary (17) and plot 
measured in £-\-A collisions as a function of the kinematic 
variables pr, yi and z. This will give a rough estimate 
of final state cff'ccts in h + A collisions. The results are 
presented in Fig. 5, and the procedure used is discussed 
below. 

Data on Rm are usually binned in either Zh, v or Q"^. 
Except for the EMC data, they are presented alongside 
the average value of the unbinned variables. For HER- 
MES data [11-13], I used the experimentally measured 
values of the DIS variables to compute the equivalent 
Pt, Ui and z. For EMC data [10], I used a computa- 
tion of the average unbinned variables from the GiBUU 
Monte Carlo generator [58, 59], which was shown to well 
reproduce the corresponding measurements at HERMES 
[60] . Another complication arises from the fact that EMC 
data have been obtained by averaging measurements at 
3 electron beam energies, Ef = 100 GeV, 200 GeV 
and 280 GeV; however, the details of such averaging are 
not immediately clear from the original paper. There- 
fore, I used the minimum and maximum of the computed 
average variables to obtain the corresponding minimum 



and maximum of the NN variables, considered as error 
band around their average value. Data have been se- 
lected to fall into 3 bins in yi (—3.0 < yi < —2.75, 
-2.5 < yi < -2.0, and -1.5 < yi < -1.1), and 1 bin in 
z {0.3 < z < 0.5). The choice of j/i bins has been made in 
order to minimize the spread of yi and z inside the bin, 
and to keep it as much as possible balanced around the 
middle value. The chosen z bin is the richest in measured 
data. Furthermore, data with similar pt and from the 
same target have been combined, with an error band in 
both Pt and Rm corresponding to the highest data plus 
error value and lowest data minus error value, the central 
value being placed in the middle. Of course, this proce- 
dure is a poor man's substitute for direct experimental 
binning in pT and yi . 

The results of Fig. 5 clearly show the evolution of final 
state cold nuclear quenching with rapidity: the quench- 
ing increases with decreasing yi . This was expected from 
the kinematic analysis of Section II, which shows a de- 
creasing v with decreasing rapidity. The size of hadron 
quenching is not small, especially for large nuclei and 
small yi rapidity. Its evolution with z is not shown in 
the figure because of large overlapping error bars aris- 
ing in the rebinning procedure. However, the original 
HERMES and EMC z^-distributions clearly show an in- 
creasing quenching with increasing Zh, especially at large 
Zh^O.5, where most of hadron production in h + A col- 
lisions takes place. Note also that quenching increases 
with the target atomic number. 

As evident from Fig. 5, the pt range covered by HER- 
MES and EMC is quite limited compared to the px for 
which hadron production in h + A and yl -I- A can be mea- 
sured. As remarked in Section III A this situation can 
be improved with a rebinning of EMC data, or with new 
measurements of hadron attenuation in ^+A collisions at 
the COMPASS experiment, which can in principle reach 
up to PT « 8 - 10 GeV. 



B. Theoretical estimate for h -\- A collisions 

As already remarked, in DIS, one has experimental 
control over all the kinematic variables. In h-|-A colli- 
sions and, most importantly for our considerations, 
z are not experimentally accessible. The non-trivial cor- 
relation of these variables with the measurable ones is 
clearly seen in Fig. 3. Moreover, the dependence of 
hadron quenching on the target atomic number A does 
not seem to follow any simple law [13, 58, 61]. For 
these reasons, it is not possible to directly use the re- 
sults of Fig. 5 to estimate cold nuclear matter effects in 
h-|-A collisions, but we need to resort to model compu- 
tations. There exist 2 classes of models which can repro- 
duce nDIS data: (i) energy loss models [61-66], which 
assume that partons hadronize well outside the target 
nucleus, and loose energy because of gluon radiation in- 
duced by rescatterings inside the target; (ii) prehadron 
absorption models [57, 58, 60, 61, 67-69], which assume 
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FIG. 5: Cold jet quenching in d+A collisions on light and heavy targets at y\ ~ —1.3, —2.25, —3.0 and z ~ 0.4, obtained from 
HERMES and EMC data on heavy and light targets. See main text for details. 



that a colorless prehadron is produced inside the target 
and can be "absorbed" via inelastic scatterings on the 
nucleons. As already remarked, the order of magnitude 
of the parton lifetime has not yet been experimentally 
or theoretically established, and both classes of models 
remain viable [18]. Hadron production in h + A colli- 
sions has a large contribution from gluon fragmentation, 
but this process has not been incorporated in absorption 
models, so far. Therefore I chose to use energy loss mod- 
els for our estimate. In particular, I will use the BDMS 
framework as implemented by Salgado and Wiedemann 
in [70, 71] and apphed to nDIS in [61, 62]. In this model, 
the nucleus is considered at rest. A parton, created with 
energy E fa v in the hard interaction, travels through the 
nucleus and experiences multiple scatterings and induced 
gluon bremsstrahlung. Hence, it starts the hadronization 
process with a reduced energy E — AE where AE is the 
energy of the radiated gluons. The reduced quark energy 
at the time of hadronization is translated into a shift 
of z in the vacuum fragmentation function D [72]. The 
medium modified FF is then computed as 



D'}^^{z, E,Eh;r}^ / dAE p{AE; Cj,, R) (27) 



1 



1 - AE/E 



AE/E 



where the quenching weight ■p(A_B) = p{AE) +pq5{AE) 
[71] is the probability distribution of an energy loss AE, 
with p(AE) its continuous part and po the probability 
of no energy loss. The quenching weight is computed for 
a static and uniform medium with characteristic gluon 
energy uJc = 0.5qL^ and size parameter R = tOcL, with 
L the medium length and q the transport coefficient of 
the medium, which characterizes the average transverse 
momentum squared gained by the parton per unit in- 
medium path- length [73, 74]. However, the nucleus den- 



sity is static but non-uniform, hence the dependence of D 
on the parton production point r, which on the r.h.s. is 
implicit in the definition of suitable static-equivalent uic 
and R [70], see Eqs. (32)- (33). They depend on a single 
parameter, the transport coefficient qo at the center of 
a reference nucleus. The outlined energy- loss model can 
well describe light hadron suppression in £ + A collisions 
at HERMES with go = 0.5 GeVVfm, fitted to tt+ pro- 
duction on Kr targets [18, 61]. I will use the same value 
for computations in ft, -I- A collisions. 

The mean free path for a parton in the target nucleus is 
A = {apAir})^^ , where a is the partonic cross-section and 
PAi'f^ the nuclear density. Assuming a to be independent 
of the atomic number, I can define a position-dependent 
transport coefficient, 



qA{b,y) = —pAib,y) 
Po 



(28) 



where qo = g^(0, 0) is the transport coefficient at the 
center of a reference nucleus of atomic number A, and 
Po — Pa(O'O)- Next, consider a parton produced at 
f — (tt, fa) which propagates in the nucleus along the r3 
direction. Its average path-length La can be defined as 



J/^ ds{s-r3)pA{rT,r3) 
dspA{rT,r3) 



(29) 



and the average nuclear density pA seen by the quark as 

ds pAirT,r3) 



PA(r) 



LA{rT-,rh) 



(30) 



Then, from Eq. (28), the average transport coefficient 
experienced by the quark can be defined as 



qA{r) = —pA{r) . 
Po 



(31) 
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For a uniform hard-sphere of nuclear density PA{r) = 
Po0{Ra — IrD^^the above definitions give La = Ra — '"3, 
Pa = Pa, and qA = as it should be. Finally, the average 
characteristic gluon energy u)c and size parameter R can 
be defined as follows: 

1 - - f°° 
^c{r) = -^QAirjL^ir) = J ds {s - y)qA{rT, s) (32) 



R{f) = a;e(r)LA(r) 



Ir^ dsqA{fT,s) 



(33) 



These equations have also been used in Rcf. [76, 77] for 
computations of jet quenching in the hot nuclear medium 
created in A+A collisions. Note that they depend on only 
one parameter, ^o- We can also see that 



2 f°° 
yy— - / ds{s-r3)qAirT,s) 



(34) 



as in Rcf. [70]. In that paper it was proven that one 
can approximate the quenching weight for a dynami- 
cally expanding medium with the quenching weight for 
an equivalent static (and uniform) medium character- 
ized by the average qA- However, the natural param- 
eters of the quenching weight arc the giuon character- 
istic energy and the size parameter. Hence, the scal- 
ing law is more properly expressed by saying that the 
equivalent static and uniform medium is characterized 
by the average Uc and R [77]. For a parton propagating 
in a static but non-uniform medium, as in our case, the 
spatial non-uniformity is equivalent to a time evolution 
of the medium. Therefore, as a rough ansatz, we may 
generalize the SW scaling law to the case of the static 
but non-uniform medium encountered in nDIS, and use 
Eqs. (32)- (33) in the quenching weight evaluation. Note, 
however, that the suitability of a single parameter q to 
describe cold nuclear matter has been recently questioned 
in Ref. [44]. 

The parton production cross-section can be computed 
in LO pQCD as discussed in Section HIB. Then the 
hadron production cross-section including cold nuclear 
jet quenching can be written as 



dp'^dyidy2 



^ dyi dy2 



(35) 



fi 



where up to terms of order 0{M'^/s), the target rest 
frame parton and hadron energy are 



E ^PT cos\i{yl) Eh = ruhT cosh(y^) , 



(36) 



where y^ = yi + log{y/s/M) are the target rest frame 
rapidities of the parton and the hadron. Isospin correc- 
tions related to the target nucleus have been included 
in the partonic cross section da^^^^^^ . Finally, we can 
quantify cold matter final state energy loss effects by the 



ratio of the above discussed cross section for collisions on 
2 targets of atomic number A and B: 



R)s{PT,y) 



daP^' 



dp7pdyidy2 



daP^~ 



■ hX 



dp'^dyidy2 



(37) 



and the amount of hadron quenching by 1 — . 

The computed R'j^ for charged pion production with 
no centrality selection is presented in Fig. 6. When ex- 
amining these plots, one should keep in mind that they 
are intended only o show the extent of the suppression 
effects on hadron production due to cold nuclear matter 
energy loss of the final state parton. They do not include 
the related transverse momentum broadening nor initial 
state effects, which will be commented on in the next 
section. The plots show a substantial final state hadron 
quenching already for midrapidity hadrons at SPS and 
FNAL energy, and for j/^j < — 2 at RHIC. The quenching 
at RHIC is reduced when increasing the rapidity, but is 
still non-negligible at yh = 0, where it is of order 5% at 
Pt > 10 GeV. This may explain the small tt" quenching 
apparent in recent midrapidity PHENIX data [41, 43]. 
Final state cold quenching at RHIC should then quickly 
disappear at forward rapidity. At the LHC, we may ex- 
pect negligible final state effects at ?//, > 3 because of 
the rapidity shift Ay in Table III. I also found a small 
hadron flavor dependence at small PhT, not shown in the 
plots, showing less suppression for kaon and proton pro- 
duction than for pion production. It would be interesting 
to compare these estimates, obtained in the quenching 
weight formalism of Salgado and Wiedemann [71], with 
the results of other energy loss implementations such as 
the twist-4 formalism of Refs. [64-66] and the reaction 
operator approach of Rcf. [44]. A nice comparison of 
the available formalisms has been recently discussed in 
Ref. [78]. 



V. INITIAL VS. FINAL STATE EFFECTS 

Before discussing the phenomenological relevance of 
the estimate of cold nuclear matter effects obtained in 
the last section, we need to discuss the importance of 
initial state effects, so far neglected. 

The initial state parton suffers multiple scatterings 
and medium-induced gluon radiation. In a simple phe- 
nomenological model [33], the resulting energy loss may 
be accounted for by a shift of the incoming parton frac- 
tional momentum, xi^xi(l — e), with e = kA^^^ the 
fractional IS energy loss. The effect of such energy loss 
is felt in a kinematic region where the flux of incoming 
partons varies rapidly with xi, typically at large rapid- 
ity. Numerical estimates from [33] indicate that IS state 
energy loss in d + Au collisions at ^ =19.4 become rel- 
evant only at forward rapidity y — ycm }t 0. According 
to the rapidity shifts listed in Table HI, we may expect 
a similar conclusion to hold for y — ycm ^ 2(5) at RHIC 
(LHC). 



11 



"T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 



0.8 



I- 



0.6 



0.4 



FNAL p+W/p+Be 

SPS p+Pb/p+p 

_l I I I I I I I I I I I I I I I L_ 



2 4 6 8 
PhT [GeV] 



0.8 



I/. 



0.6 



0.4 



I I I I — I I I I I I I I I I I I I I I . 



/ 



/ RHIC d+Au/p+p 

y.=o J 

yh=-i 

y,=-2 

yh=-3 

I I I I I I I I I I I I I I I I I I I 



5 10 15 20 
PhT [GeV] 



FIG. 6: Energy loss model estimate of final state hadron quenching in cold nuclear matter for midrapidity pions at SPS and 
FNAL, and several negative rapidities at RHIC. 



If the final state parton is long lived, as assumed in the 
theoretical estimates of the previous section, the medium 
affects hadron production mainly through elastic and ra- 
diative energy losses. In this case, the FS energy loss 
enters the computations as a shift in z of the fragmen- 
tation function, see Eq. (27). Hence, differently from IS 
energy loss, it is large in regions where the fragmenta- 
tion functions changes rapidly in z, namely at large z. 
At fixed p/iT the average {z) increases with decreasing 
rapidity and decreasing (see Fig. 7 right). Coupling 
this with a decrease in final state parton energy Ej with 
decreasing rapidity, it is easy to explain the behavior and 
large size of final state suppression shown in Fig. 6. 

A consistent framework for considering the interplay of 
IS and FS energy loss in the reaction operator formalism 
is discussed in Ref . [44] , which presents numerical results 
for the partonic fractional energy loss AE/E in a case 
study of a homogeneous medium of fixed length L ^ 5 fm. 
At any given parton energy E, the FS fractional energy 
loss is generally smaller than the IS fractional energy loss. 
They both start at around 10% when £; = 10 GeV, but 
FS energy loss tends rapidly to as -B increases, contrary 
to IS energy loss which stabilizes around 5% at E > 1 
TeV. However, particle production at fixed rapidity in h+ 
A and ^-|-v4 collisions shows a strong correlation between 
the IS parton energy Ei and the FS parton energy Ej, 
see Fig. 7 left. As a result, for midrapidty hadrons at 
SPS we have comparable IS and FS state energy loss of 
5-10%. For midrapidity hadrons at RHIC, FS energy 
loss becomes quite small, and IS radiation is about 5%. 
For backward rapidity production, yh — Ucm = —3, FS 
energy loss is now larger than IS energy loss, viz., 10% vs. 
5%. A detailed computation including realistic nuclear 
geometry is needed to quantify their effect on hadron 
spectra. 

In summary, IS and FS cold nuclear matter effects are 
expected to be dominant in different rapidity regions. 



viz., at forward and backward rapidity, where the es- 
timates presented in this paper and in Ref. [33] indicate 
that they are large. Their effect on the midrapidity region 
has to be more carefully and quantitatively considered: 
it depends on the center of mass energy of the collision, 
and can be expected to decrease with increasing y^. 



VI. DISCUSSION AND CONCLUSIONS 

In this paper, I have examined the role of final state 
interactions in cold nuclear matter in modifying hadron 
production on nuclear targets with leptonic or hadronic 
beams. Initial state parton energy loss has been consid- 
ered in [22, 23, 33]. Since in £ + A collisions only FS 
interactions are present, I built a kinematic dictionary 
that relates the variables used for the discussion of nDIS 
and heavy-ion A + B collisions, and demonstrated the 
(limited) extent to which available experimental data on 
hadron suppression in nDIS can give direct information 
on final state cold nuclear matter effects in A + B colli- 
sions. In this respect, the EIC [50, 51] will be able to efh- 
ciently cover the regions in phase space which the HER- 
MES [13] and EMC [10] experiments could not examine. 
A nearly full coverage of the SPS and FNAL phase space 
may alternatively be achieved either by a reanalysis of 
EMC data, or by new measurements of hadron attenu- 
ation at the COMPASS experiment [52]. The latter op- 
tion is particularly interesting: COMPASS has a similar 
kinematic coverage to EMC, but higher luminosity and 
very good particle identification capabilities. Therefore 
a, fj. + A program at COMPASS, building on the knowl- 
edge accumulated at the HERMES [13] and CLAS ex- 
periments [14, 79], would greatly improve our knowledge 
of the space-time evolution of hadronization, and gather 
vital data for the interpretation oi h + A and A + A col- 
lisions and the quest for the Quark Gluon Plasma. 
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FIG. 7: Correlation of initial and final state parton energy (left), and hadron transverse momentum and fractional momentum 
(right). 



Hadron production \vl h ^ A and A + A collisions is 
affected by cold nuclear matter in 2 ways. 

1. IS and FS energy loss, and possibly FS prehadron 
absorption, suppress hadron spectra by non negligi- 
ble amounts at forward [33] and backward rapidity, 
respectively. 

2. The transverse momentum broadening associated 
with induced radiation and multiple scatterings 
in the medium will modify the hadron phT spec- 
trum, further suppressing it at ^ 1 — 2 GeV 
and enhancing it at intermediate momenta up to 
PT«5-6 GeV [27]. 

I used an energy loss model based on the EDMS for- 
malism and tuned to i + A data, to estimate the size of 
final state cold hadron quenching in hadronic collisions, 
which was found to be large at midrapidity at fixed tar- 
get SPS and FNAL energy, and at backward rapidity at 
RHIC energy. It will be interesting to compare this result 
with estimates based on the GLV [44] and high-twist [64- 
66] formalism for energy loss, and on nuclear absorption 
models [57, 58, 67]. 

In this paper, only the induced FS suppression has 
been discussed and estimated. However, for phenomeno- 
logical applications to h + A and A + A collisions the in- 
terplay of suppression and enhancement of hadron spec- 
tra will need to be elucidated. In h + A collisions, cold 
nuclear effects will be important to understand the evo- 
lution with rapidity of the nuclear modification factors, 
which decreases at forward rapidity [15-17] but increases 



at negative rapidity [28-31]. In A collisions, the time 
scale for the formation and equilibration of the Quark- 
Gluon Plasma is much larger than the time scale for cold 
nuclear matter interactions: tcoid ^ ^eg- Hence, both 
cold and hot quenching should be considered for a correct 
interpretation of experimental data. At RHIC midrapid- 
ity, I estimated cold quenching of hadron spectra to be of 
order 5-10%. It is much smaller than the observed factor 
4-5 hot medium quenching observed in central Au + Au 
collisions, and is negligible in first instance. At LHC cold 
nuclear matter effects are likely to be negligible in a large 
midrapidity interval |y| < 3. At SPS energy, where cold 
nuclear matter effects may be of the same order of mag- 
nitude or larger than hot medium effects, they both need 
to be taken into account in any QCD tomographic anal- 
ysis to detect and extract the properties of the Quark 
Gluon Plasma. 
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